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We present a measurement of the top-quark mass in pp collisions at = 1-96 TeV which uses events 
with an inclusive signature of missing transverse energy and jets. The event selection is sensitive 
to tt —> W^bW'b —> lubqq'b independent of the lepton flavor, and results in a large acceptance for 
W TV decays. AU-hadronic ti decays and events with identified electrons or muons are vetoed 
to provide a statistically independent sample with respect to all previous measurements. The top- 
quark mass is inferred from the distribution of the scalar sum of all jet transverse energies and 
the missing transverse energy. Using 311 pb^^ of integrated luminosity recorded by the Collider 
Detector at Fermilab, wc measure a top-quark mass rat = 172.3lg''g*(stat.)±10.8(syst.) GcY/(?. 
While the uncertainty on mt is larger than that of other measurements, the result is statistically 
uncorrelated with those of other methods, and thus can help to reduce the overall mt uncertainty 
when combined with other existing measurements. 

PACS numbers: 12.15.Ff, 13.85.Ni, 14.65.Ha 
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The top-quark mass, rrit, is an important free parame- 
ter in the Standard Model (SM) of particle physics. Being 
roughly 40 times larger than the mass of its weak isospin 
partner, the b quark, rrit gives large contributions to 
electroweak radiative corrections which, when connected 
to precision electroweak measurements, can be used to 
derive constraints on the masses of the yet-unobserved 
Higgs boson [ij, and of particles belonging to some SM 
extensions At the Tevatron pp collider top quarks 
are produced mainly in pairs through quark-antiquark 
annihilation and gluon-gluon fusion processes. Because 
the Cabibbo-Kobayashi-Maskawa matrix element Vtb 0| 
is close to unity, the SM top quark decays to a boson 
and a b quark almost 100% of the time. The final state of 
a top-quark pair thus includes two W bosons and two b- 
quark jets. When only one W decays leptonically, the tt 
event typically contains a charged lepton, missing_trans- 
verse energy (I^x) from the undetected neutrino a-nd 
four high-transverse-energy jets, two of which originate 
from b quarks. 

Recently the CDF collaboration has reported precision 
mt measurements using tt events containing identified 
high-pT leptons (e, /x) j3] and all-hadronic decays @. In 
this paper we describe a top-quark mass measurement 
which uses events collected by a multijet trigger, and 
selected by requiring an inclusive high-p^ neutrino sig- 
nature, consisting of large Events containing identi- 
fied high-pT electrons or muons {E^ > 20 GeV, > 20 
GeV/c), as defined in f?!, are removed in order to increase 
the relative contribution of W tv decays and provide 
a statistically independent sample with respect to other 
lepton-based measurements [Hj. All-hadronic ti decays 
are discarded by the I^rp requirement so that orthogo- 
nality with respect to the all-hadronic mass sample is 
ensured Unlike previous analyses based on the iden- 
tification of ^ ei'{fii') and W —f qq' decays, our event 
selection does not permit a full kinematical reconstruc- 
tion of the tt final state. For this reason, the top-quark 
mass is derived from the Hx distribution, where Ht is 
defined as the scalar sum of all jet transverse energies 
and the I^j^. 

Results reported in this paper are obtained using 311 
pb~^ of integrated luminosity from pp collisions at ^/s — 
1.96 TeV, recorded by the Collider Detector at Fermi- 
lab (CDF II). The CDF II detector is described in detail 
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elsewhere [8| . It consists of a magnetic spectrometer sur- 
rounded by a calorimeter and muon system. The mo- 
menta of charged particles are measured up to a pseudo- 
rapidity of I77I = 1.0 in a cylindrical drift chamber, which 
is inside a 1.4 T superconducting solenoidal magnet. Sili- 
con microstrip vertex detectors, located immediately out- 
side the beampipe, provide precise track reconstruction 
useful for vertexing and extend the coverage of the track- 
ing system up to I77I = 2.0. Electromagnetic and hadronic 
sampling calorimeters, arranged in a projective-tower ge- 
ometry, surround the tracking systems and measure the 
energy and direction of electrons, photons, and jets in the 
range \r]\ < 3.6. In addition, the good hermeticity pro- 
vided by the calorimeter allows the detection of high-py 
neutrinos by the measurement of the l^rp. Muon systems 
outside the calorimeters allow the reconstruction of track 
segments for penetrating particles. The beam luminosity 
is determined using gas Cherenkov counters surrounding 
the beam pipe, which measure the average number of 
inelastic pp collisions per bunch crossing. 

The data sample used in this analysis is collected by a 
multijet trigger which requires four or more Et > 15 GeV 
clusters of contiguous calorimeter towers, and a scalar 
sum of transverse energy clustered in the calorimeter of 
J2 > 125 GeV. The initial data sample consists of 4.2 
million events and is further reduced offline by the appli- 
cation of kinematical and topological requirements aimed 
at optimizing the tt signal significance [9;] . Briefly, we re- 
quire at least four jets having Et > 15 GeV and \ri\ < 2.0; 
l^rp significance, ifx^ , greater than 4.0 GeV^/^, where 
is defined as / -^t ; and a minimum separa- 
tion in azimuthal angle between the l^j- and the closest 
jet, min A(/) > 0.4 rad. In our selection, jets are identified 
as groups of calorimeter tower energy deposits within a 
cone of radius AR = ^/A(jP + Ar]^ < 0.4, and their ener- 
gies are corrected for calorimeter non-linearity, losses in 
the gaps between towers, multiple interactions, and par- 
ticle response calibrations [l3| ■ This selection reduces the 
data sample to 597 events, with a signal to background 
ratio S/B ~ 1/5. In order to further increase the ex- 
pected S/B ratio and reject background events with only 
light quark or gluon jets, &-quark jets ("& tags") are iden- 
tified by the reconstruction of secondary decay vertices 
using the SECVTX algorithm, as in 0]. After these selec- 
tions and the requirement of at least one 5-tagged jet, we 
observe 106 events with S/B ~ 1; about 44% of the signal 
acceptance is accounted for by tt — > W^bW~b rvbqq'b 
decays, while the remaining ti content is dominated by 
e(/i)-|-jets events, in which the lepton fails the standard 
high-pT identification cuts. 

Background events with b tags arise from QCD heavy 
flavor production, electroweak production of W bosons 
associated with heavy flavor jets, and from false iden- 
tification by the SECVTX algorithm. The overall num- 
ber of background b tags in the final data sample, and 
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FIG. 1: Observed and expected Ht background distributions 
in data control samples depleted of signal contamination, see 
text for details. 



their corresponding kinematical distributions, are esti- 
mated using a per-jet parameterization of the 6-tagging 
probabiUty derived from the multijet sample. For the 
parameterization, we use events with exactly three jets, 
having Et > 15 GeV and |?7| < 1.0, where the tt con- 
tent is negligible. The parameterization exploits the b- 
tag rate dependencies on the jet Et, the charged track 
multiplicity inside the jet cone, and the projection of the 
I^rp along the jet direction in the transverse plane, which 
is defined by 1^^'' — cos A(/)(^y, jet). The extrap- 
olation of the 3-jet 6-tagging probability to higher jet 
multiplicity events, and the capability of the parameter- 
ization to track sample composition changes introduced 
by the kinematical selection, are checked using > 4-jet 
data samples depleted of signal content, as described 
elsewhere (a) data before the optimized kinematical 
selection on and min A0(^2n, jets); (b) Ijf^ < 3.0 
GeVi/2, minA</)(^j,,jets) > 0.3 rad, and (c) > 3.0 



min A0(^2"i jsts) < 0.3 rad. As a result, the 
6-tag rate parameterization is found to predict the num- 
ber of background b tags, and the shape of their cor- 
responding kinematical distributions, to within 10% in 
the 4 < A'jet < 6 region, where 96.4% of the tt signal is 
expected after the optimized kinematical selection. Fig- 
ure [T] shows the comparison between expected and ob- 
served background Ht distributions in the data control 
samples (a) , (b) , and (c) . The expected Ht distributions 
are derived from the 6-tag rate parameterization applied 
to each jet belonging to a given data sample, before 6-jet 
identification requirements. The observed Ht distribu- 
tions receive one entry per 6-tagged jet for a proper nor- 
malization with the expectation. The normalization and 
shape of the observed and expected distributions are in 
good agreement for all control samples. 

The final data sample, after the optimized kinemati- 
cal selection and the additional requirement of at least 
one 6-tagged jet, contains a total of 127 6-tagged jets. 
The number of 5-tagged jets yielded by background pro- 
cesses in that sample is expected to be nl^^ = 57.4 ± 8.1. 
The excess in the number of b tags is ascribed to top- 
quark pair production. We derive a measurement of the 
top-quark mass from the observed Ht distribution. The 
Ht distribution from the selected data is fit to the sum 
of signal and background Ht contribution parameteri- 
zations using an unbinned likelihood technique. Proba- 
bility density functions are determined for signal, as a 
function of rrit, and for background events by fitting a 
functional form from the corresponding Ht distributions 
(templates). For consistency with our per-jet background 
prediction method, the Ht distributions from data and 
simulated signal events receive one entry per 5-tagged jet. 

We calibrate our method using events with inclusive tt 
decays generated with different input values of rrit rang- 
ing from 150 to 200 GeV/c^, in steps of 2.5 GeV/c^. 
These events are simulated using the herwig [ii| gener- 
ator in conjunction with the CTEQ5l [T^j parton distri- 
bution functions (PDFs), QQ 13] for the modeling of b 
and c hadron decays, and a full simulation of the CDF II 
detector [l^. They are then subjected to the same 
selection as the recorded events. The Ht distributions, 
derived at discrete values of the top-quark mass, are pa- 
rameterized by a continuous functional form as a func- 
tion of rrit in order to smooth the distributions and in- 
terpolate between the templates. For any given rrit the 
probability to observe a particular Ht value is specified 
by a normalized Pearson type IV function , in which 
the parameters are assumed to be linearly dependent on 
rrit- The best parameterization is determined by a si- 
multaneous binned likelihood fit to all signal templates. 
In Fig. [21 four signal templates are shown overlaid with 
their corresponding parameterization. 

For background, the Ht distribution is derived from 
the 6-tag rate parameterization applied to jets belonging 
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FIG. 3: The background Ht template, after the subtraction 
of the content (using mt=172.5 GeV/c^), is shown overlaid 
with the fitted parameterization. 



FIG. 2: Four Ht signal templates with mt ranging from 150 
to 195 GeV /c? . Overlaid are the fitted parameterizations at 
each generated top-quark mass. 



to the kinematically selected data sample, before &-jet 
identification requirements. It has no dependence on the 
top-quark mass, except from a procedure adopted to sub- 
tract the expected signal content (~ 15% for mt = 172.5 
GeV/c^). The arbitrary rat choice in the subtraction 
procedure is accounted for in the background shape sys- 
tematic uncertainty. A single probabihty density func- 
tion, defined as the sum of a gamma function and two 
gaussians, is used to fit the background Ht template, as 
shown in Fig. [3l 

The likelihood function used to extract the top-quark 
mass includes as free parameters the number of expected 
signal and background h tags (rij and n?,), and mt- It is 
specified by three factors: 



bkgi 



where 



N 



C-sh{mt) = Jl 



Us + Hb 



ris+rib 



^bkg — 



iV! 



N 



(1) 



(2) 

(3) 
(4) 



and N is the number of observed b tags in the final data 
sample. In Cshijnt) the product is over the number of 



observed b tags, and Psig{H^\mt) and Pbkg(-ffT) are the 
probability density functions for signal and background, 
respectively. The second factor of Eq. ((T|) represents a 
Poisson constraint on the total number of b tags observed 
in the data. Finally, in Eq. ^ the background normal- 
ization is constrained to its expected value n^'^^ to within 
(T„^ = 10% • n1^^ ■ The likelihood is maximized with re- 
spect to 71s, nb and mt- The statistical uncertainty from 
the fit procedure is taken from the mt values where the 
log-likelihood changes by 0.5 units from its maximum. 
Since we are counting b tags and not events, the Ht dis- 
tribution is not strictly Poisson in nature. We correct for 
this effect below. 

We use simulated data ensembles (pseudo- 
experiments) to check our fitting procedure for possible 
systematic biases. For each generated top-quark mass 
from 150 GeV/c^ to 200 GeV/c^, we construct pseudo- 
experiments, with the same statistical properties as our 
observed data sample, by randomly sampling from the 
signal and background templates. Then we perform like- 
lihood fits to each pseudo-experiment and characterize 
the accuracy of the technique in determining the correct 
mt value. In each pseudo-experiment, the number of 
background b tags is Poisson fluctuated around its 
expectation, n1^^, while the number of signal b tags 
is Poisson fluctuated around the number observed in 
the data, minus the central value for the background 
expectation. In this procedure, b tags from single and 
double 6-tagged events are fluctuated separately. For 
each pseudo-experiment, the likelihood flt provides 
the measured mt along with the positive and negative 
statistical uncertainties from which pull distributions 
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(GeV) 



TABLE I: Relevant sources of systematic uncertainty. 



Source 


Amt (GeV/c^) 


Jet energy scale 


9.6 


Generator 


3.8 


Background shape 


2.1 


PDFs 


1.5 


ISR 


0.9 


FSR 


0.9 


Background fraction 


0.8 


6-jet energy scale 


0.7 


Trigger efficiency 


0.7 


Limited Monte Carlo statistics 


0.6 


b tagging 


0.5 


Total 


10.8 



FIG. 4: Ht distribution from the selected data sample, over- 
laid with the expected distribution from the unbinned likeli- 
hood fit. The inset shows the —2\nC from the final fit as a 
function of rat- 



are derived. The mean of the puU distribution, averaged 
as a function of the input mt, is consistent with zero 
(—0.01 ± 0.02), while the width is slightly larger than 
unity, due to the inclusion of duplicated Ht values 
in the pseudo-experiment distributions in the case of 
double-tagged events. For the current analysis, we 
correct for this effect by scaling the statistical errors 
taken from Aln£ = —1/2. The scale factor is the pull 
width averaged over mt ranging between 150 and 200 
GeV/c2, giving 1.08 ± 0.02. 

Applying our method to the observed Ht distribu- 
tion, we find Ug = 76.2 ± 11.4, n?, = 54.6 ± 5.1, and 
mt = 172.3t9°g®(stat.) GeV/c^. The statistical un- 
certainties on mt are consistent with expectation from 
pseudo-experiments performed with an input top-quark 
mass of 172.5 GeV/c^. The result from the fit to the data 
is shown in Fig.[31 The inset shows the function — 21n£ 
from the final fit as a function of mt ■ 

Systematic uncertainties arise from uncertainties in our 
understanding of the detector response and in the as- 
sumptions employed to infer the top-quark mass from 
the observed data. For each source of systematic uncer- 
tainty, the relevant input quantities are varied by ztlcr, 
and new signal or background Ht templates are pro- 
duced by performing the event selection and reconstruc- 
tion on the modified samples. Then these new fixed tem- 
plates are used to run pseudo-experiments. The mean 
shift in the fitted top-quark mass with respect to the 
input value is taken as the systematic uncertainty associ- 
ated with the given assumption or effect. Table U reports 
all the relevant sources of systematics associated with 



our measurement. The dominant source of uncertainty 
(9.6 GeV/c^) given the choice of Ht as discriminant vari- 
able is associated to jet energy scale (JES) uncertainty. 
For each jet considered in the Ht calculation the relative 
jet energy scale uncertainty, which is mainly driven by 
uncertainties on particle response calibrations and the 
out-of-cone jet energy modeling, varies from 3% to 8% 
depending on rj and pt of the jet. We determine the 
impact of the JES uncertainty on our measurement us- 
ing pseudo-experiments in which the nominal jet ener- 
gies are varied by ±1 standard deviations. Additionally, 
the dependence on the Monte Carlo generator is esti- 
mated as the difference in the extracted top-quark mass 
in PYTHIA [l3| and HERWIG events, and amounts to 3.8 



GeV/c^. Other sources of uncertainty are related to the 
background shape and normalization and are evaluated 
to be 2.1 and 0.8 GeV/c^, respectively. We estimate the 
uncertainty from PDFs using signal samples in which the 
events are weighted according to their probability to oc- 
cur using different sets of PDF eigenvectors. The system- 
atic uncertainty is computed by considering differences 
between the CTEQ5l and mrst72 pJi] PDFs parameter- 
izations, different Aqcd values, and the sum in quadra- 
ture of half the difference between the ±lcr shift of the 20 
CTEQGm uncertainties, for a total of 1.5 GeV/c^. Varia- 
tion of initial (isr) and final state (fsr) gluon radiation 
settings, as in are found to contribute 0.9 GeV/c^ 
of systematic uncertainty each. Systematic uncertainties 
due to the 6-jet energy scale, trigger simulation effects, 
statistically limited Monte Garlo samples, and 6-tagging 
efficiency modeling, are small and give a combined er- 
ror of 1.2 GeV/c^. The total systematic uncertainty is 
estimated to be 10.8 GeV/c^ assuming all sources to be 
uncorrelated. 

In conclusion, we report the first top-quark mass mea- 
surement using inclusively selected ^j^-l-jets tt events 
with a large acceptance for W ^ tv decays. The re- 
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suit, mt = 172.3tg°6^(stat.)±10.8(syst.) GeV/c^, is com- 
plementary and statistically independent with respect to 
precision CDF measurements 0| , and consequently, al- 
though not competitive by itself, it will help to reduce by 
a few percent the overall uncertainty on mt when com- 
bined with other existing results. 
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